Magnetic field controlled vacuum charge in graphene quantum dots with a mass gap 
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The effect of a magnetic field on the charged vacuum is investigated. The field dependence of 
the energy levels causes jumps in the total vacuum charge that occur whenever an energy level 
crosses the Fermi level and this leads to re-entrant cycles of vacuum charging and discharging. In 
atomic systems these effects require astrophysical magnetic fields of around 10 s T but in graphene 
with a mass gap they occur in laboratory fields of about 1 T or lower. It is suggested that an 
electrostatic graphene quantum dot defined by a gate electrode provides a solid state model of the 
as yet unobserved charged vacuum as well as a model of an atomic system in an extreme astrophysical 
environment. Phase diagrams are computed to show how the total vacuum charge depends on the 
confining potential strength and applied magnetic field. In addition the field dependence of the 
vacuum charge density is investigated and experimental consequences are discussed. 

PACS numbers: 73.22.Pr, 73.21.La, 31.30.J- 



The long-sought-for charged vacuum [J Q is the 
ground state of strong field quantum electrodynamics 
(QED). Usually the vacuum is neutral but it charges in 
the presence of an electric field strong enough to lower 
a bound state into the negative energy continuum. For 
example, when the charge on the nucleus of a hydrogenic 
atom increases to beyond ~ 172 the Is state enters the 
negative energy continuum and if this happens at con- 
stant Fermi level, the vacuum charges. This is accom- 
panied by spontaneous emission of two positrons which 
would enable the effect to be observed if a strong enough 
field could be created. In principle, this is possible be- 
cause the critical charge can be exceeded in a collision 
between two uranium nuclei. But the interaction time is 
too short to allow the transition to a charged vacuum to 
occur [3] and, despite much effort, vacuum charging has 
not yet been observed. However, it may be possible to 
observe it in a semiconductor analogue. 

The band gap of a semiconductor is analogous to the 
mass gap of an atomic system and semiconductor quan- 
tum dots are analogous to natural atoms In electro- 
static quantum dots, electrons are confined by an electro- 
static potential that is generated by a gate electrode and 
replaces the Coulomb potential of a natural atom. Nor- 
mally the dot is engineered so that the confined electron 
energies are just below the edge of the conduction band. 
But in a material with a small band gap it should be 
possible to use a stronger potential to lower a state into 
the valence band and create a charged vacuum analogue. 

Any material that allows a state to be lowered into the 
valence band with a modest gate voltage is suitable. One 
candidate is monolayer graphene on substrates that in- 
duce a gap, for example BN Q, Ru [f| and controversially 
SiC ■ Other candidates include semiconducting carbon 
nanotubes and narrow gap semiconductors. Graphene is 
the only candidate with a Dirac-like energy dispersion 
and hence the candidate that provides the most precise 



analogue of QED. Indeed it has already been suggested 
that a charged vacuum occurs in graphene in thepres- 
ence of a Coulomb impurity with enough charge [8| but 
it is difficult to vary the impurity charge experimentally. 
However any attractive potential of sufficient strength 
charges the vacuum. This means a graphene quantum 
dot is an accurate and practical analogue of the atomic 
charged vacuum. In addition, dots made of materials 
with non-relativistic energy dispersion may allow studies 
of unusual charged vacua whose properties differ from the 
atomic one. 

But the most interesting feature of the quantum dot 
analogue is that the charged vacuum is extremely sensi- 
tive to an external magnetic field. A rough estimate of 
the field needed to cause significant effects may be ob- 
tained by equating the rest mass-energy to the cyclotron 
energy. For atomic electrons this gives about 10 10 T, an 
ultra strong magnetic field that only occurs in extreme 
astrophysical environments such as magnetized neutron 
stars P. In contrast, for graphene with m c 2 ~ 100 meV 
and c ~ 10 6 ms _1 the same estimate gives about 10 T, 
well inside the laboratory regime. 

Moreover, important effects already occur at lower 
fields. Depending on their quantum numbers, energy lev- 
els both rise and fall as the magnetic field increases. If the 
energy of a charge carrying state rises above the Fermi 
level, the vacuum discharges while it charges if the en- 
ergy of a state falls below the Fermi level. This leads 
to re-entrant cycles of vacuum charging and discharg- 
ing. The relevant energy scale for these processes is the 
depth of the state in the continuum. For a hydrogenic 
atom with Z — 172 this is about 15 keV [2| and the field 
needed to discharge the vacuum is about 10 8 T, still in 
the astrophysical regime, while for a graphene quantum 
dot it is < 1 T. These processes do not seem to have 
been investigated before and are studied here in the con- 
text of a graphene quantum dot. But the wave equation 
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for a graphene dot is just the Dirac equation so exactly 
the same physics should occur in atomic systems in ultra 
strong magnetic fields. The graphene dot is not only an 
analogue of the charged vacuum but also an analogue of 
atomic physics in extreme astrophysical environments. 

The objectives of this Letter are, first to demonstrate 
that charged vacuum states occur in graphene quantum 
dots, secondly to demonstrate magnetic field induced vac- 
uum charging and discharging and finally to consider the 
experimental consequences. The dot is taken to be cir- 
cularly symmetric and the electrostatic potential is mod- 
elled by V(r) = Vo exp(— (r/A) p /2) where r is the radial 
co-ordinate, |Vb| is the well depth and A = 50 nm is its 
width, p determines the shape and slope of the well; the 
bottom flattens and the edge sharpens has p increases 
(Fig. [TJ top left). The magnetic field, B, is taken to 
be uniform and perpendicular to the dot plane and the 
quantum states are found by solving the two dimensional 
effective mass equation. Interestingly, the reduced spa- 
tial dimensions may make it easier to realise a charged 
vacuum The mass is generated by including a site- 
dependent splitting parameter in the Hamiltonian, the 
same approach as in earlier work on graphene with a 
mass gap d, [ll[ . 

The effective mass Hamiltonian for graphene con- 
sists of two 2x2 blocks which together are equiva- 
lent to the 4-component Dirac Hamiltonian. One block 
gives the states near the K point of the Brillouin zone 
and the other gives the states near K' . The states 
near K are obtained from the 2-component Hamiltonian 
H = (7/?i)er • (p + eA) + V + m,oc 2 <j z , where the a 
are Pauli matrices, p is the momentum and A is the 
magnetic vector potential. Here c = j/h, 7 is taken 
to be 646 meV nm [l2j and mnc 2 is taken to be 100 
meV, the upper end of the observed range (10 - 100 
meV). The eigenstates of a circularly symmetric dot are 
4>(r) = (xi( r ) exp(i(m— 1)0), X2M exp(im0)), where 9 is 
the azimuthal angle and m is the total angular momen- 
tum quantum number. Equations for the radial functions 
are obtained by making the substitutions f± — \frx\ 1 
if 2 = y/rX2- This leads to 
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These equations are discretized with a second order 
forward-backward difference scheme [l3| which leads to 
a real, symmetric eigenvalue problem that is solved nu- 
merically. The states near K' are found in a similar way. 
The numerical system radius is 600 nm. 

Physically, the valley index is equivalent to a pseudo- 
spin and the states near K and K' correspond to pseudo- 
spin up and down respectively. The total angular mo- 
mentum h(m — 1/2) is the sum of the orbital angular 
momentum and pseudo-spin. The orbital angular mo- 
mentum quantum number I is therefore to — 1 for states 
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FIG. 1: Model potentials (top left); well depth and field de- 
pendent energy levels (bottom) and typical states at 0.704 T 
(top right). 



near K and to for states near K' . 

Figure [TJ shows the typical behaviour of energy levels 
and quantum states. The lower left frame shows the K 
energy levels as a function of well depth when B = 0, 
m = 1 (i.e. I = 0) and p = 2. The results clearly show 
bound state levels plunging into the negative energy con- 
tinuum as the well depth increases and similar behaviour 
is found for other values to and p. The lower right frame 
shows the K energy levels as a function of B for the same 
parameters as in the lower left frame and well depth 265 
meV. The energy levels in the negative energy continuum 
[lij move down as B increases but those of the bound 
states that have entered the continuum move up. The 
bound state level closest to the continuum edge crosses 
E = —moc 2 (dashed line) when B ~ 2.4 T and this cor- 
responds to vacuum discharge in a system with Fermi 
level, Ef ~ — moc 2 . The top right frame shows typi- 
cal states. When a bound state enters the continuum it 
hybridizes with the continuum states and forms a Fano 
resonance [lj|. Individual states that contribute to the 
resonance closest to the continuum edge are shown in 
the figure. The resonance width depends on the strength 
of the hybridization. Semi-classical analysis [l6| shows 
that a forbidden region surrounds the dot. As shown by 
the arrows in Fig. [TJ the width of this forbidden region 
decreases when p = 8 and this strengthens the hybridiza- 
tion [l7|. When p — 2, the resonance width is less than 
the numerical continuum level separation (~ 1 meV) and 
the resonance consists of one state. But when p = 8, it 
involves about 2-4 states. In the case of Fig. [JJ B is cho- 
sen so that /i has the roughly the same amplitude for 
the 2 main contributing states and the resonance width 
is about 3 to 4 meV. The superposition of bound and con- 
tinuum character is clearly visible. In contrast, strongly 
hybridized states do not normally occur in Coulomb po- 
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FIG. 2: Total vacuum charge as a function of well depth (left) 
and magnetic field (right). 
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FIG. 3: Vacuum charge phase diagrams for p = 2 (left) and 
p = 8 (right). K phase boundaries solid, K' dashed. 



tentials because the width of the forbidden region is large. 

The vacuum charge density is the charge density 
induced in response to an external potential, V. It may 
be found from the commutator QED charge operator, 
(— e/2)[ip, 7 C V], where 7 is a Dirac matrix, ip is the 
Dirac field operator and ip is its adjoint. Or it may be 
found from the normally ordered QED charge operator, 
N(—e^pj°ip). Alternatively it is given by the standard 
form p(V) — p(0) where p is found by summing over states 
below Ef- The two QED operators are identical [l8j . 
Further, both are equivalent to the standard form and 
this follows from their expectation values. The vacuum 
expectation value of the commutator operator is p(r) = 
(-e/2)(£ £ „ <BF ,J</w(r)| 2 - £ B „>^,J<^w(r)| 2 ), 
where a represents the component and valley indices. 
The standard form is -eJ2 E <e f a l^«"( r )l 2 ~ P(°) = 

(-emZ En <E F , a \<t>n a (r)\ 2 - E En >E F , a \<t>n a (r)\ 2 ) + 

[(-e/2)(E £ „< £F ,J^(r)| 2 + E En >E F , a \<f> : n a (r)\ 2 ) - 
p(0)]). The terms in the square brackets vanish because 
of completeness and chiral symmetry [l9j so the stan- 
dard form is equivalent to the QED forms, p is used to 
find the charge density and its integral gives the total 
vacuum charge Q = (-e/2)(£ lEn<Ep<a -T,E n >E P<a )- 
In QED the sums in p and Q are divergent and have 
to be treated with charge renormalization but they are 
convergent in the present numerical model because the 
energy spectrum is bounded. 

The total vacuum charge as a function of well depth 
and magnetic field is shown in Fig. [2j Real spin split- 
ting is included and the effective g-factor is taken to be 
2.0. E E is just above — m c 2 so all real-spin-split Landau 
levels remain below it. For B — the vacuum charge in- 
creases monotonically with well depth in a series of steps. 
The first step occurs when the I = level shown in Fig. [1] 
enters the vacuum. Each step has height 4 and this cor- 
responds to a 2-fold pseudo-spin degeneracy and a 2-fold 
real spin degeneracy. At constant well depth the vacuum 
charge as a function of B shows re-entrant behaviour for 
both values of p. For instance, for well depth 220 meV 
and p = 2, it falls when B ~ 0.41 T and then rises when 
B ~ 3.9 T. The fall occurs because the I = 0, K level 
enters the continuum while the rise occurs because the 
I = — 1, K ' level leaves the continuum. The behaviour at 
larger well depth is similar but richer because some levels 



go through a minimum as a function of B [19j. Another 
effect of the field is real spin splitting. This allows an 
odd numbered vacuum charge and leads to double steps 
of height 1, for example at B ~ 3.9 T with well depth 
220 meV and p = 2. 

To investigate the details, vacuum charge phase dia- 
grams are computed (Fig. [3]). Each phase boundary in- 
dicates where an energy level, E n crosses the Fermi level 
and is given implicitly by E n (Vo,B) — E E . The total 
vacuum charge is shown on selected portions of the di- 
agrams. For clarity real spin splitting is not included. 
Hence each line corresponds to a vacuum charge step of 
2 electrons when B ^ and 4 when B = 0. The re- 
sults in Fig. [2] are sections through the phase diagrams 
with real spin splitting included. The phase boundaries 
reflect the physics of the system: they have have small 
splittings at B = 0, unless I — 0, and there is pronounced 
B-dependent splitting. These effects can be understood 
from the non-relativistic limit of the effective mass equa- 
tion. To order 1/toq the splitting at B — results from 
the pseudo-spin-orbit interaction and the independent 
splitting results from the interaction of the pseudo-spin 
with the magnetic field. Quantitatively, the exact split- 
tings at B = Q, V = -180 meV and p = 2 are 1.05, 1.88 
and 2.49 meV for the lowest states at I = 1, 2, 3 respec- 
tively, while the lowest order pseudo-spin-orbit splittings 
are 1.36, 2.57 and 3.63 meV. The -B-dependent splitting 
at 1 T, calculated from the pseudo-spin (/-factor is 6.33 
meV while the exact splitting at I = is 6.05 meV. Hence 
the non-relativistic limit describes the physics qualita- 
tively but the exact equation is needed to find the split- 
tings accurately. 

Vacuum charge densities are shown in Fig. [3J The 
charge density increase associated with the charging steps 
is found from p(Vo+SVo) — p{Vo) where SVo is small (—0.1 
meV). States below the threshold for vacuum charging 
are described as undercritical while those above it are 
called overcritical [2[ . The left side of the figure shows the 
charge density associated with the first step at B = T 
in Fig. [21 To a good approximation, the vacuum charge is 
stored in one overcritical state when p = 2 and two when 
p = 8. The open circles indicate the charge density com- 
puted from these few states. Near the peak they agree 
with the exact data to better than a few parts in 1000. 
The exact density is also well approximated by the under- 
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FIG. 4: Comparison of vacuum charge density and over and 
undercritical state densities (left); magnetic field dependence 
of p = 2 vacuum charge density for 4 and 8 electrons (right). 

critical bound state density (crosses). The ith overcriti- 
cal state may be approximated [2| by + 5Z c i>icV'c> 

where ipb is the undercritical bound state and tp c are un- 
dercritical negative continuum states. The coefficients an, 
and bi C form a unitary transformation Q hence the ap- 
proximate vacuum charge density reduces to \ipb\ 2 - The 
numerical data shows this approximation is accurate to 
a few parts in 1000 near the peak. A consequence of this 
approximate sum rule is that the vacuum charge den- 
sity is hardly affected by hybridization with continuum 
states, although individual states contributing to it are. 

The right hand frames of Fig. @] illustrate the effect of 
the magnetic field on the vacuum charge density when 
p = 2 and the total charge is 4 and 8 electrons. A 
magnetic field normally compresses the charge density 
because the cyclotron length is proportional to l/y/B. 
However, in the present case the density expands because 
states of higher orbital angular momentum enter the vac- 
uum when the B increases. These states have a larger 
spatial extent so the density expands. For example, with 
4 electrons at T the charge density is composed of the 
4-fold degenerate I = state while at 5 T it is composed 
of real and pseudo spin split states with / = and I = — 1. 
The expansion continues until all the higher angular mo- 
mentum states are exhausted. 

There are several experimental consequences of these 
findings. First, the vacuum charge in a graphene quan- 
tum dot may be detectable via emssion of holes that 



is analogous to spontaneous positron emission. It may 
be possible to pump the gate voltage to enhance this 
effect. Secondly, the strong hybridization may be de- 
tectable with scanning tunnelling spectroscopy (STS). 
Although the hybridization does not affect the density, 
the states in Fig. Q] extend to a radius about an order of 
magnitude larger than the dot radius. They should be 
detectable when the experimental resolution is less than 
the width of the Fano resonance and this requires tem- 
peratures of a few K. Thirdly, it may be possible to detect 
the vacuum charge directly for example by capacitance 
measurements or quantum point contact electrometers. 
The ideal material for these experiments should be un- 
doped or p-doped to ensure that an empty state is taken 
through the gap and it should be graphene with a gap to 
ensure that the effective mass equation is a precise ana- 
logue of the Dirac equation. Existing results |f| H[ sug- 
gest that a suitable material can be found. In addition 
it should be possible to observe a charged vacuum ana- 
logue in narrow gap semiconductors or semiconducting 
carbon nanotubes, although the effective mass equation 
is no longer relativistic. 

In summary, a quantum dot in graphene with a mass 
gap provides an accurate model of the charged vacuum. 
The vacuum charge is strongly affected by a magnetic 
field and a sufficiently high magnetic field can discharge 
the vacuum. These effects should also occur for atomic 
electrons in magnetic field but the small rest mass-energy 
of graphene charge carriers, allows them to be seen with 
fields ~ 1 T instead of the extreme astrophysical fields 
needed for atomic electrons. Experimentally, it may be 
possible to detect the graphene charged vacuum directly, 
or by observation of hole emission or by probing the states 
with STS. The effect of interactions remains to be inves- 
tigated but a Thomas-Fermi approach @ suggests the 
atomic vacuum charge survives the effect of interactions. 
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